Modeling Pionic Fusion 

Alexander Volya, Scott Pratt and Vladimir Zelevinsky 
National Superconducting Cyclotron Laboratory, Michigan State University, 
East Lansing, Michigan 48824-1321, USA 
(February 9, 2008) 



Recently observed rare heavy ion fusion processes, where the entire avail- 
able energy is carried away by a single pion, is an example of extreme collectiv- 
ity in nuclear reactions. We calculate the cross section in the approximation 
of sudden overlap, modeling the initial and final nuclei by moving harmonic 
oscillator potentials. This allows for a fully quantum-mechanical treatment, 
exact conservation of linear and angular momenta and fulfillment of the Pauli 
principle. The results are in satisfactory agreement with data. Mass number 
CN , dependence and general trends of the process are discussed. 
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Nuclear fusion reactions which produce a pion are often referred to as pionic fusion. 
Pion production has been observed [1|-|3| at energies approaching absolute threshold, where 
the entire available energy is converted into the pion, demonstrating an amazing collective 
behavior of nucleon systems. However it remains quite difficult to incorporate the observed 
collectivity into existing theoretical models. A variety of studies [§-0] have dealt with 
subthreshold pion production in heavy ion collisions, where the energy per nucleon is below 
the energy threshold of the elementary single-particle reaction NN — ► NN+tt. Most models, 
such as those featuring pion bremsstrahlung mechanisms f|[|, quantum molecular dynamics 
approaches ||10|| , perturbative calculations using Boltzmann-Nordheim-Vlasov equations (Tl| 
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or nuclear structure functions ||12|| , provide a good picture at energies starting from E/A « 30 
MeV up to the single-particle threshold E/A = 280 MeV in the laboratory frame. In the 
present work motivated by the experimental results of |l|-|3| our aim is to consider even 
lower energies and study the behavior of the cross section of fusion reaction in the region 
down to ~10 MeV above the absolute threshold. This necessitates a careful consideration 
of limitations on the reaction given by conservation laws and the Pauli exclusion principle 
which govern the behavior of the cross section in this extreme situation. The statistical 
approach used in most existing models at higher energy has to be substituted by low energy 
many-body structure physics. 

Our model, that is described in the Sect. [TT], considers the cross section in the Born 
approximation, assuming that pion production occurs through coupling to a single nucleon. 
All possible further rescatterings of the pion are expected to significantly reduce the proba- 
bility of the reaction, and are ignored as higher order processes. A schematic picture of the 
reaction is shown in Fig. [I] demonstrating the pionic fusion of two nuclei A and A'. 
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FIG. 1. The pionic fusion of two nuclei in the sudden approximation is illustrated. 



Many-body nuclear mean field parameters are assumed to be constant and suddenly 
change from the initial to the final values. This will be referred to as the sudden ap- 
proximation. The three-dimensional harmonic oscillator shell model is used to describe 
the structure of the incoming and outgoing nuclei. This allows analytical calculation of 
all necessary overlap amplitudes. The stationary wave functions are constructed as Slater 
determinants projected onto good angular momentum. Taking into account the center-of- 
mass motion we preserve linear momentum. Sect. |TJ shows the implementation of the 
model for the case of pionic fusion of two identical nuclei. In Sect. [TV] we present a low pion 
momentum approximation, for which more general results could be derived. The parallel 
discussion of mathematical details is given in the Appendices. The application of this model 
to experimentally observed pionic fusion reactions shows a good agreement with data. The 
comparison is presented in Sect. |V| along with some predictions for heavy nuclei. 



II. DESCRIPTION OF THE MODEL 



A. The transition amplitude 



We approach the problem of pionic fusion as a stationary scattering problem. We consider 
the reaction cross section to be given by the Fermi golden rule in terms of the transition 
amplitude (F\H\I), where I and F refer to the initial and final states, respectively, of 
the whole system including the emitted pion. The density of final pion states is given by 
Vk 2 dk dQ/(2ir) 3 with fcasa momentum of the pion produced; Q is a solid angle in the center 
of mass (CM) frame, and V stands for the quantization volume. In all further calculations 
the pion is assumed to be fully relativistic whereas nucleons obey non-relativistic quantum 
mechanics. We use a set of natural units with h = c = 1. In this framework the differential 
cross section can be written as 
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where m is a nucleoli mass, p n is the CM momentum per nucleon in the initial state and 



to = yk 2 + m\ is pion energy. 

On the single-nucleon level one can use a phenomenological Hamiltonian density for the 



pion-nucleon interaction [|i~3|j, 



7i = gt/j'fsTipjT + An — — ipn ■ nip + An — ■ n x nip . (2) 
m n m* 

A number of studies have been performed analyzing this form of the interaction within the 
context of chiral perturbation theory ||14|| . The first term, often called in the literature the 



impulse or Born term, is responsible for single-pion production in a p-wave. We neglect the 
second and the third s-wave terms which require an additional interaction to absorb the extra 
pion created in the four-point vertex. We believe that due to the difficulty of recombining 
the nucleons into an appropriate final state the second and third terms become increasingly 
unimportant for larger nuclei. It has also been experimentally observed that in the pionic 
fusion reactions the pion is predominantly produced in the p-wave . Reduction of the 
first term in the Hamiltonian to a non-relativistic case gives an interaction of the form 

a • k 

with the coupling g appropriately defined according to isospin. Separation of the quantized 
pion field, 

n(x) = £ ^^(4e-* k - x + a k e + ^\ (4) 

in the matrix element of Eq. reduces the transition amplitude to the following form 

<^> = 7^^-</l E 9k-ae~^} , (5) 
V2uV 2m nuclcons 

where |/) and \i) are final and initial states of the nucleon system. 



B. Nuclear wave functions 

We will approximate a state of a nuclear system with an antisymmetric combination built 
upon single-particle (s.p.) states. We take these states from the harmonic oscillator shell 
model, which allows for the analytic calculation of corresponding overlaps. The approach 
however can be extended to any single-particle basis. Each of the single-particle states 
can be characterized by the number of excitation quanta in three Cartesian directions, the 
nucleon spin and isospin projections. The locations of the centers of the harmonic oscillator 
potentials for all separate nuclei have to be introduced as additional parameters to the wave 
function. The importance of these parameters in projecting a nucleon wave function onto a 
state with correct total momentum for every nucleus participating in the process is discussed 
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below. Following these assumptions we will write the wave function of a nucleon system as 
follows 

(A+l)th s.p. state 

c * * 

I Si, £1; a 2 , s 2 , t 2 ;...;r), (a A+i , s A+1 , t A+1 ; a A+2 , s A+2 , t A+2 ; r )) . (6) 

S v ' S v ' 

nucleus A nucleus A' 

In this example we assume that the system consists of two nuclei A and A' with the centers 
of their respective harmonic oscillator potentials at r and r'. The single-particle orbitals are 
numbered from 1 to A for the first nucleus and from A + 1 up to the total number of nucleons 
Af = A + A' for the second one. Labels ) are Cartesian quantum numbers 

of single-particle states, while s and t are the spin and isospin projections. Protons and 
neutrons can be considered separately as well as different spin projections of the nucleons, 
reducing the wave function of the state to a product of four components. If the described 
separation is performed and the resulting part of the wave function contains only single- 
particle states with the same values of either s or t then the corresponding index is omitted in 
writing. We use a standard form for the one-dimensional harmonic oscillator wave functions 
centered at r in the coordinate representation: 



(x\(a;r) v ) = ^-^H a (v(x-r))e-^- r ) 2 / 2 . (7) 

The parameter v is defined for a single oscillator as v = ^Jmu. These parameters characterize 
the mean field potentials for every incoming or outgoing nucleus. The function H n (x) is the 
nth order Hermite polynomial of the variable x. The discussion of the overlap integrals such 
as ((a' ,r') v i\(a,r) v ) , and the general form of the results is presented in Appendix A. 

A simple projecting technique was used to construct wave functions as eigenstates of the 
momentum operators that correspond to the total momenta of each individual nucleus, 

|(ai, s x , p), (d A+ i, s A+1 ,t A+1 ; p')) = 

AT 1 / [ + °°d 3 rd 3 r' \(a 1 ,s 1 ,t 1 ;.. .;r),(a A+1) %i ) tA + i;...;r')>e i(p ' r+pV) . (8) 

J J -co 

It is easy to check that 

A 

Vj|(ai,...;p),(aA+i,- ■•',p')) = P • • • ! p), («a+i, • • • ; p')) ( 9 ) 

i=i 

and 

A+A' 

~ i Z! Vj|(ai,...;p),(aA+i, ...;p')) = p'l(«i,- --IP), ((Sa+i, • • ■ ; p')> ■ ( 10 ) 

j=A+l 

In the above example the situation with two-nuclei state is shown, which is appropriate for 
describing the initial state in pionic fusion. The final state containing just one fused nucleus 
is constructed analogously. 

Due to the finite range of the interaction, the overall normalization M of the state @ 
that contains several moving nuclei, is just a product of normalizations for each of the 
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constituent nuclei individually. It is useful to write the CM coordinates separately from the 
relative coordinates of the nucleons 



|(ai; d 2 ; .. .d A ;r) v ) = \(d C u = (0,0,0);^^) • (11) 

The relative coordinate wave function \iprei) can be complicated but the CM part for the 
ground state nucleus is simply represented by the unphysical ground state oscillation of the 
center of mass in the effective harmonic potential with the parameter v\f~A. This is removed 
by a projection (|§|) onto the correct momentum state. The normalization integral can be 
expressed as 



M 2 = j J d 3 rd 3 r'((d C M;r') vVJ \(dcM;v) vVX ) ^ re # re i)e ip(r - r,) 

= J d 3 rj rf 3 rV A( r - r ')^/Vp(r-r0 = f^j^y^/A^ {u) 

A different method of calculating the normalization along with further justification of this 
form for the CM part of the wave function is discussed in Appendix B. We also note here 
that with a slight modification of Eq. ( |T!?D the orthogonality of the nucleon wave functions 
can be shown 

((di, d 2 ; . . . d A ;p')v\(di, d 2 ; . . .d A ; p)„) = A/" 2 5 P; p/ . 



III. FUSION REACTIONS A + A -> 2A + vr 

For the remainder of the paper we will assume A to be the mass number of each of the 
initial nuclei with proton-neutron composition (Z,N) and w the oscillator parameter. The 
entire initial state is characterized by a set of the single-particle quantum numbers {di}. 
The fusion product has 2 A = Af nucleons, the oscillator parameter v, and the final state 
quantum numbers {A}. The collision is considered in the CM reference frame; therefore 
we use p and — p to denote the momenta of the incoming nuclei and k for a final pion 
momentum with corresponding P/ = — k as the total momentum of the recoil nucleus. The 
integration of the wave functions leading to correct momenta, Eq. (^), is performed at a 
final stage so initially overlaps are calculated as functions of r , r' and R, the locations of 
the centers of the two initial nuclei and the final nucleus, respectively. 



A. Charged pion production 

We begin with the case of tt + production where one of the initial protons interacts with 
the pion field producing a neutron and a real on-shell pion. With the assumption that 
the pion was produced in a single-particle interaction, the total amplitude of the process 
becomes a sum over all possible amplitudes shown in Fig. ^|, with the pion vertex connecting 
any of the initial protons to any of the final state neutrons with the correct relative sign to 
preserve antisymmetry. 
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FIG. 2. One of the amplitudes of the total fusion process: an initial proton from the nth orbit 
produces a tt + and ends at the Ith. final neutron single-particle orbit. F n is the remaining overlap 
of a proton system with the nth initial single-particle state missing. G[ is the neutron overlap with 
no Ith state in the final system. 

Suppose the interacting proton in the single-particle state n produced a neutron in the 
state I of the final nucleus. In the initial state we sum over the occupied orbitals of the first 
and the second colliding nucleus, for n < Z and for Z < n < 2Z, respectively. We use the 
notations Gi(r,r' , R) for the neutron overlap 

Gi = §i, . . . /9j_i,s,_i;/3, + i,sj +1 ; . . . ; R) |(ai,si; . . .;r)(3 A +i, sa+i] ■ ■ • ;r')) , (13) 

S v ' 

no Ith s.p. state 

F n (r, r', R) for the proton overlap 

F n = ((A, si; . . .; R)| (a 1; s x ; . . . ; a n _i,s n _i; a n+1 , s n+ i; ...)), (14) 

V v ' 

no nth s.p. state 

and H n i for a single-particle matrix element 

f ((A,S /; R)|(7a-k e - ik ' x |(a n , Sn ;r)) n<Z 

^n^ = <^ . • (15) 

1 ((A, S z ; R)|p^ • ke- ik - x |(«„, s n ; r')) n > Z 
Finally, following Eq.(KI), the total amplitude can be expressed in terms of the following sum: 



(F\H\I) = k } Kr J_ — / / d 3 rd 3 r'd 3 R V (-l) n+i F n G /J ff nZ e ip{r - r ' ) - ip ^ R . (16) 



The determinants of the matrices are constructed from a product of the single-particle 
overlaps of size (2N) x (2N) for the neutrons (Gi) and (2Z — 1) x (2Z—1) for the protons (-F n ). 
The Gaussian nature of the single-particle overlaps allows one to separate all exponential 
factors that govern the general trend of the cross section leaving only some polynomials of 
a general form that carry spin, isospin and Pauli blocking information. These mathematical 
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manipulations are discussed in some detail in Appendix B. Here we present a final expression 
for the square of the transition amplitude 



K*W)| 2 = %(^\ l2 ^^ A - l) |M + | 2 , (17) 



V 2 \v 2 Aj urn 2 

in which the exponential factor £, effective oscillator parameter r\ and reduced amplitude 
M + are introduced as follows 

/ 2p 2 k 2 k 2 (w 2 -v 2 ) \ 

4-exp^ ^ 2 V 2 + W 2 2Av 2 (v 2 + w 2 )) ' ^ ^ 

V=^~ - , (19 

and 

M+ = P(fc, p)e p - k/Al ' 2 + Q{k, p)e~ p - ]l/Av2 . (20) 

Here, P and Q are dimensionless polynomials of p and k, the total CM momentum of the 
initial nuclei and the final pion momentum. The polynomials are to be determined using 
particular configurations of the initial and final nuclei. They are also functions of v and w 
which determine the appropriate momentum scale. If the two colliding nuclei have the same 
initial shell model state then P(k,p) = ±Q(k, —p) (the phase difference given by ± sign for 
even or odd Z, respectively, is due to imposed Pauli antisymmetry, see Appendix B). The 
procedure of analytically calculating P and Q involves finding the determinant of the matrix 
constructed from polynomials that result from integrating a product of Hermite polynomials 
of the form ((/3, R)\(a, r)); and performing the integrational Fourier-type conversion, Eqs. 
(H). This process is discussed in Appendix B. The size of the matrices is determined by the 
number of nucleons of the same spin-isospin type. 



B. Neutral pion production 

The case of 7r° production can be considered in a similar fashion. A neutral pion can 
be produced either by one of the protons or by one of the neutrons, which couple with a 
negative relative sign. Compared to charged pions the coupling is larger by a factor y/2. 
The final amplitude can then be expressed, similarly to Eq. (0), as 

l< F ^ J >l 2 = £ (t^) 3/2 l M °l 2 • ( 21 ) 

V 2 \Av z J lorn 1 

Here the reduced amplitude can be split into a proton and a neutron part: 

M = P p (k,p)e p ^ Av2 + Q p (k,p)e- p ^ Av2 - P n (k,p)e p ^ Av2 - Q n (k,p) e - p ^ Av2 . (22) 
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IV. LOW PION MOMENTUM APPROXIMATION 



Due to the specific form of the polynomials discussed above, further simplifications can 
be made for the case of ir° production. Near the absolute threshold, the pion momentum 
|k| is small compared to all other momentum parameters |p|, v and w, and can be ignored 
in polynomials. Then 



H n i = (((3i, si;TL) v \a ■ ke 



- ik-x I 




« ( Sl \a ■ k\s n ) ((A; R)„|(a w ; v) w ) exp ^ - ^ - j . (23) 

With this approximation, the interaction part is factorized into exponents as shown in the 
expression above. Therefore the total pionic fusion amplitude is a product of a pure fusion 
amplitude and the expression that arises from the operator a ■ k acting on the nucleons. 
For a given type of the initial and final nucleon, the sum of a single-particle matrix element 
multiplied by the corresponding overlap of the remaining particles reduces to a sum of matrix 
elements multiplied by the corresponding minor which is related to a determinant of a full 
matrix. It is shown in Appendix C that the polynomials can be expressed in an analytical 
form if all inner harmonic oscillator shells are completely filled without any gaps in all 
participating nuclei. This restriction allows any type of particle-hole excitations within the 
outer unfilled shell. 

The total differential cross section for a neutral pion production close to absolute thresh- 
old is given in the form: 



da = g 2 Ak 3 (2n_\ 3/2 6{A „ 1)+Qf+Qi (4w\ Q ^ 2pW 
dVl (2Tr) 2 2pm \Av 2 J 1 \AvJ 




2 1 ~ 2 
- 2 M . (24) 
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Here the integers qj, j = x,y,z are introduced as differences between numbers of quanta 
in final and initial systems for three Cartesian directions; Qi and Qf are total numbers of 
quanta in initial and final systems. These values are defined as 

Qj = Pj- Yj a i' Q = H (a x + a y + a z ), Q f = {Px + Pv + P*)- 

nucleons nucleons nucleons nucleons 

(25) 

The spin and radial parts of the wave function are completely decoupled in our non- 
relativistic description of the nucleon system. This allows to introduce the matrix element 
used in Eq. (||) 

^=]H</I £ 7-,^-k|5>, (26) 



nucleons 



where i and / are the spin-isospin parts of nucleon wave function of initial and final systems, 
respectively. This matrix element could be directly computed for every particular nuclear 
configuration, but for a large number of states degenerate within harmonic oscillator model 
it is useful to use an approximation for the average 
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M = (Z^-Zt-Ni + Nt). (27) 

The Cartesian directions of the harmonic oscillator quantization axes are chosen in such a 
way that the z axis coincides with the beam direction, though the spin is quantized along 
the k axis that simplifies the action of a ■ k which is used to obtain Eq. fl2"7|). Integers Z^, 
Z[, N<i and iVj are mean numbers of particles for each spin-isospin combination with respect 
to our axis of spin quantization. The polynomials T n (x), defined in Eq. fl36|) of Appendix 
A, can be approximated as 

^ W33 - W£T ■ (28) 

This approximation is valid in the limit that the arguments become large and allows for 
a better quantitative understanding of the behavior of the cross section. The value of the 
argument is almost independent of the mass number A at threshold energy: 

P\ -7— « 6. 

V Ar/vw 

In Eq. ( p4|) only the lowest order term in the pion momentum is retained resulting in 
a p-wave cross section (exponents with k are also ignored). The equation has only one 
numerical parameter 7, the origin of which is discussed in Appendix C. This parameter is a 
product of four factors, one for every spin/isospin nucleon species. Each factor depends on 
the number of particles of corresponding type and on their distribution within the highest 
harmonic oscillator shell for both initial and final nuclei. Numerically, 7 range from 1 to 10 
for light nuclei. The cross section can be zero if some symmetries are not preserved (spin, 
isospin, oscillator symmetry) as well as by virtue of Eq. ( |39l) in Appendix A if creation of 
the final system requires an odd number of quanta relative to the initial system in any of 
the transverse directions. 



V. APPLICATION OF THE MODEL AND RESULTS 



A. The reaction p + p — ► d + tt + 

The first and the simplest example to calculate is the two-nucleon fusion reaction p+p — > 
d + 7r + . This example serves here only for illustrative purpose as we do not include pion 
rescattering due to the full interaction given by Hamiltonian density of Eq. (0) which is 
important for this elementary process. Moreover, the deuteron hardly can be approximated 
with the harmonic oscillator shell model. The polynomials P and Q in this case do not 
depend on p being equal to the matrix element of a ■ k/|k| evaluated between the spinors 
of initial interacting proton and final neutron. In Eq. (|T7D we choose a minus sign for 
antisymmetry. Here, P and Q correspond to the choice of the first or second initial proton 
to produce a pion, respectively. 

Dominant partial wave channels are summarized in the following table along with our 
results for their reduced amplitudes. The table was constructed by separation of initial 
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singlet and triplet states of the NN system. Partial waves of the ir — d system printed in 
the left column that yield the dominant contributions to the amplitudes which are shown in 
the right column. 



pion NN state amplitude 



s-wave 3 P 1 2 V2 smh{^f± 
p— wave 1 5'o 



p— wave 



(29) 



2 cosh(^f^; 



As can be seen from the table above, this cross section is predominantly p-wave in nature 
at low pion energies. The s-wave contribution that comes |Tj| from rescattering of the pion 
due to the interaction (0) was not included. The total cross section averaged over spin 
projections in the initial state and summed over final states is 



da g 2 k 3 

d£l 2mp\f2ixv 3 



exp 



Ap 2 + k 2 
2v 2 



3 cosh ( 2k P™ sd ) _ ! 



The obtained p-wave cross section behaves at low energies as 

a(pp — > dir + ) = a (k/rrij,) 3 , 



where 



a 



2\/2-ng I m 



5/2 



-2mm n /v 2 



rn 



3/2y2 



(30) 



(31) 



(32) 



Choosing the oscillator parameter v = 216 MeV/c reproduces the experimental value |15 
45" m 0.42 fm 2 . For this case the fusion is sensitive to the tail of the wave function in 
momentum space. Since the wave function of a deuteron is extremely non-Gaussian with a 
long tail in coordinate space, choosing v to reproduce the deuteron's r.m.s. charge radius 
would result in a grossly underpredicted cross section. For the fusion of heavier ions, the 
incoming nuclei are moving at a slower velocity and their momenta per nucleon are similar 
to characteristic momentum scales of the wave functions. 

The oscillator parameter v can be best obtained by matching used here harmonic oscil- 
lator type deuteron wave function to its experimentally known behavior Hl~5| . The choice of 



this parameter between 180 and 220 fm would lead to the values of 4a in the range of 0.06 
to 0.48 fm 2 . 



B. The reaction 3 He + 3 He -> 6 Li + 7r+ 

As a next step, we apply the model to the experimentally studied pionic fusion reaction 
3 He + 3 He — > 6 Li + tt + , where even first excited states of the 6 Li nucleus have been resolved 
0. This reaction involves heavier nuclei so that the process of pion rescattering becomes less 
important as discussed above. The polynomials P and Q for Eq. flT7| ) can be constructed in 
a direct way considering the shell model structure of all nuclei involved in the reaction. The 
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ground 1 + and first excited 3 + states of 6 Li were constructed within the P3/2 j-subshell. In 
Fig. the total cross section for this reaction is calculated for the fusion into the ground 
state (left panel) and the first excited state (right panel). The contributions of the s-wave 
and p-wave to the cross section are plotted together. We choose a value v = 118.91 MeV/c 
for 6 Li as it corresponds to the oscillator frequency of 15.06 MeV, the parameter of MK3W 
model The initial parameter w = 112.7 MeV/c was chosen by assuming the r.m.s. 

size 2.14 fm of 3 He. In Fig. |4] we show the differential cross section for this fusion reaction 
going into the ground state of 6 Li (solid line) and the first excited state (dashed line). The 
beam energy is assumed to be fixed so that the corresponding absolute values of the pion 
momentum are 96 and 90 MeV/c, respectively. 

Comparison with the experiment || in which pionic fusion resolves the few lowest levels of 
6 Li shows that we obtain a reasonable ratio of the cross sections. However, we underpredict 
the magnitude by approximately 40%, compared to the estimated experimental value of 
111 ±11 nb for the ground state transition. We note that the result is sensitive to parameters 
of the shell model and their choice in the harmonic oscillator approximation is quite uncertain 
for light nuclei. For example, a variation of the final oscillator frequency within 10% range 
of the used value would lead to the values of the cross section between about 20 and 140 nb. 
Using more realistic non-Gaussian wave functions might significantly improve the model. We 
might also be underestimating the cross section due to inherent limitations of the approach. 
For instance, we do not consider a gradual change of the nuclear mean field in the process 
of fusion substituting it with the sudden approximation. 




50.0 0.0 

k„ (MeV) 



50.0 100.0 



FIG. 3. Reaction cross sections for 3 He + 3 He — ► 6 Li + ir + . The left panel shows the transition 
to the ground state and the right panel to the first excited state of 6 Li at 2.18 MeV. The solid lines 
represent the total cross section, dashed and dotted lines are s and p-waves, respectively. 
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FIG. 4. Differential cross section of the reaction 3 He + 3 He —* 6 Li + ir + . On the left panel the 
solid line represents the transition to the ground state of 6 Li and the dashed line to the first excited 
state; the corresponding absolute values of pion momentum are 96 and 90 MeV/c, respectively. 
The right panel displays the experimentally observed values [2] of the differential cross section of 
the transition to the ground state (squares) and to the first excited state (circles) of 6 Li. 



C. The reaction 12 C + 12 C 24 Mg + tt' 



Here we apply our approach to the cross section of the 12 C + 12 C — > 24 Mg + tt° reaction. 
This process, along with its isospin analog 12 C + 12 C — > 24 Na + tt + , represents those few 
heavy ion pionic fusion reactions for which experimental data exist [0. The application 
of the developed theory does not present a great difficulty except the fact that the cross 
section is quite dependent on the structure of initial and final states of interacting nuclei. 
Within the harmonic oscillator picture we have approximately 3 x 10 8 different combinations 
of interacting states that correspond to the same energy. Angular momentum and isospin 
conservation constraints reduce this number by several orders of magnitude. Additional 
shell model interactions have to be introduced to build up a realistic nuclear state for each 
of the nuclei and reduce this large number of states, that are degenerate in our model, to 
the ones of interest. Based on this argument we will present here the Monte-Carlo averaged 
cross section, where we average over random Cartesian states. In the following Fig. [| we 
display the total reaction cross section as a function of pion momentum. We use here the 
oscillator parameters v = 104 MeV/c and w = 119 MeV/c which are estimated by various 
theoretical models [|T?|j[0I • 
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FIG. 5. The reaction cross section for C + C — > Mg + tt with oscillator parameters 
v = 104 MeV/c and w = 119 MeV/c as a function of pion momentum. 

The experimentally estimated cross section for this reaction is 208 ± 38 pb which was 
observed for pion momentum 41 MeV/c |T]|. In this example we again underestimate the 
cross section. To see the sensitivity of our results we present in Fig. || the dependence of 
the cross section on oscillator parameters for pion energy at about 6 MeV (momentum 41 
MeV/c). 




FIG. 6. The total cross section of C + C — > Mg + tt as a function of the model parameters 
v and w. The calculation is done for a pion momentum 41 MeV/c which corresponds to the total 
energy of about 6 MeV above threshold. 

This figure indicates that a reasonable variation of parameters could cause a change in 
the answer by an order of magnitude. We emphasize again that in our calculations we did 
not project the participating nuclei onto appropriate shell-model states. Such a projection 
would require additional nuclear structure input. Given that the existing experimental data 
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do not clearly resolve the structure of the final state this seems sufficient. As a conclusion, 
within all the limitations discussed above, the agreement between the introduced theory and 
the experimental results of this rare process seems to be remarkable. 



D. Calculations for heavy nuclei 

In this section we apply the low-momentum approximation for the cross section described 
by Eq. ( f24f) to several reactions, with the goal of understanding the general dependence 
with respect to the mass of the incoming nuclei. In order to calculate the cross section, one 
needs the harmonic oscillator parameter v which can be estimated from the experimentally 
determined r.m.s. radii of the nuclei [ 20(1 , 



1 ^ / ,\ 1^1/ 3 



7 E<>f> - 3 E£(« + i)- (33) 

In order to calculate the cross section, one needs to know the incoming energy of the 
nuclei as well as the energy of the outgoing pion. Calculations of the cross sections were 
performed for incoming nuclei 9 Be, 12 C, 16 and 20 Ne with corresponding fusion products 
18 0, 24 Mg, 32 S and 40 Ca in the limit of low pion momentum. In this limit the cross section 
is proportional to the cube of the pion momentum, 

a = a(k 3 /ml). (34) 

Values of a are displayed as a function of the mass number of the incoming nuclei in Figure [7]. 
The shell model configurations are again randomly chosen from the available set of Cartesian 
states that conserve isospin and parity. Average values are represented by filled diamonds 
while the states with the highest and lowest cross sections are represented by the boundaries 
of the error bars. The large error bars demonstrate the wide fluctuation in strengths for 
individual states. However, despite the fluctuations, it is clear that the overall trend is of a 
decreasing cross section with increasing mass. 

Also shown in Figure [7| are experimental measurements represented by open circles for 
the pp, 3 He 3 He and 12 C 12 C cases discussed previously. The corresponding calculations, 
which were performed for the experimentally measured pion momenta rather than in the 
low-momentum limit are also displayed with closed circles. One sees that the cross sections 
fall by several orders of magnitude, but the measurements are still feasible throughout the 
wide range of masses. Calculations could be performed for heavier nuclei, but for larger 
masses the Coulomb barrier becomes important, and shuts off the possibility of fusion for 
masses larger than 20. 
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FIG. 7. The general behavior of the pionic fusion cross section A + A — > 2A + n versus the mass 
number of initial nucleus A. The plotted value a is related to a total cross section as a = a {k/m^f". 
Calculations in the low-momentum limit (filled diamonds) show that cross sections fall by several 
orders of magnitude in this mass range, but remain in the picobarn region for nuclei as large as 
oxygen. The highest and the lowest cross section found within the shell model configurations are 
represented by error bars. Experimental measurements are displayed (open circles) and compared 
to calculations (filled circles) which were performed for the finite pion momenta corresponding to 
the experiments. 



Near threshold meson production represents a unique area of heavy ion reactions. In 
this area the reactions underline the pronounced features of quantum many-body physics. 
Most theoretical approaches to understanding and predicting these phenomena lose their 
validity in such an extreme regime. In this paper we have proposed a simple model to 
study the processes of deep subthreshold pion production. The pionic fusion cross section 
was obtained in a Born approximation with respect to pion production and in the sudden 
approximation for the nuclear rearrangement. The participating nuclei were described by 
the harmonic oscillator shell model in moving oscillator potentials. The advantage of the 
method is that it allows one to incorporate energy, momentum, spin and isospin conservation 
laws precisely and respect the Pauli principle at all steps of the calculation. Further aspects 
of nuclear structure could be additionally taken into account. At threshold energies these 
constraints pose the most powerful restriction on the reaction and cannot be ignored as is 
done in statistical and kinetic models which are reasonable at higher energy. The obvious 
disadvantage of the model is that the sudden approximation does not consider the slow 
changes of the nuclear mean field in the process of interaction. For the future developments 



VI. CONCLUSIONS 
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it seems feasible to incorporate the time dependence and solve the equations for the evolution 
of a nuclear mean field parameters. 

The nearly analytical form of overlaps greatly simplified the calculations for this study. 
We used a spherically symmetric nuclear mean field but in some cases this symmetry pro- 
hibits the transition and this would require a consideration of deformations, i.e. different 
oscillator parameters for different directions. The above mentioned limitations are reflected 
by the difficulty in determining the parameters of the model, and lead to about an order 
of magnitude ambiguity in the result for s-d shell nuclei. More realistic single-particle wave 
functions could be incorporated into the model. Some of the exponential factors in Eqs. ( |TTD 
and (|2ll) arise directly from the Fourier transformation of the Gaussian tails in harmonic 
oscillator wave functions and could be substituted with modified factors that would reflect 
a more realistic behavior. 

We would like to stress here again that pionic fusion is a very rare process presenting 
a tiny fraction of the total cross section. The agreement that was observed between cal- 
culations and experimental data for the cross sections ranging from 10~ 4 to 10~ 9 barns is 
remarkable. Within the limits of the low pion momentum approximation in the class of the 
reactions A + A — > 2 A + tt, we were able to obtain a general formula, Eq. (p4|) , for the cross 
sections. The proposed techniques can certainly be applied in the same manner to other 
pion production reactions. The processes of electrofission |21] present an another interesting 
avenue to exercise this technique. 
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APPENDIX 



A. Harmonic oscillator wave functions and overlaps 

Work with harmonic oscillator wave functions often involves integrations of the expres- 
sions constructed of polynomials and Gaussian weights. Thus, the following integral is useful 

(a > 0), 

[ +0 °x n e- ax2+ibx dx= x [^e~ b2 / ia a - n l 2 T n {^\ , (35) 
Joo V a \2^/a J 

where T n (x) is a sum arising from the binomial expansion, 

' (36) 

This expression can be used for the evaluation of any integral encountered in this work. 
There are two important limiting cases for the sum T n (x), x — > and i>1: 
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(n-1)!! ■( 

v 9n/ 2 it n is even 

\imT n (x) = \ , Jim T n (x) = 2~ n/2 x n . (37) 

if n is odd 



The Gaussian- Fourier integration of Eq. (|35| ) is a transformation on the space of polynomials 
defined on the basis 

x n - [* n ](p) = r„(p) (38) 

The following two-dimensional integrals often appear in our calculations, 

na~ l ~ n l 2 {n — 1)!! if n is even 
(x-y) n e- a{x2+y2) dxdy= { , (39) 

if n is odd 

/ J{x-y) n e-< x2+ ^e l ^ x -^dxdy=^{^ ll \- p2 l 2a T n . (40) 



The basic block of the calculations is the overlap of two one-dimensional harmonic os- 
cillator wave functions with different oscillator parameters, shifted locations of the centers 
and possible additional factor e~ tkx that enters the single-particle interaction integral from 
Eq.(||). This type of integral, the generalized Debye- Waller factor, can be written in a 
factorized form: 

mr) v \e-* kx \(a ] r') w )=r ] 1 / 2 x 

c::p I ~ fc2 (r-r')W ik{rv 2 + r'w 2 ) \ {{r ,_ r)k . vw) m 

eP \2(v 2 + w 2 ) 2 (v 2 + w 2 ) (v 2 + w 2 ) ) M[T r ^' fc ' u '^' ^ 

where r\ is given by Eq. (|i~9| ) and P nm (r, k;v,w) is a dimensionless polynomial of r and k of 
the highest power n + m with coefficients dependent on w and v. The following are examples 
of these polynomials for the smallest values of n and m: 



Poo(r, fc; v,w) = l 



P 01 (r,k;v,w) = P 10 (-r, k;w,v) 



,2> 



y/2 (ik + 

V 2 + W 2 

2vw {y 2 + w 2 — (k — ir v 2 ) (k + irw 2 )) 



Pn(r,k;v,w) — — — . 

The technique of obtaining these expressions is simple though tedious. An important situa- 
tion k = would correspond to the overlap of two wave functions without a pion production, 
in this case we will not write k as an argument. It can be shown that |I9[ 



k =^ti fiH^ (-l) j ~ l v i - l w i - k (rn\ i+j ~ k ~ l 
P^(r;v,w)= E \ V^- n, 4) P k i(0;v, W ) (42) 



fc+/=0,2,4.. 

and 



k\ i\ (i-k)\ (j-iy. Vv^, 



Pki(0;v,w) = JjVli- + f / 2 2 (v) , (43) 



17 



with Pf being the associated Legendre polynomials. A bit simpler case is 

I -| -| min(i,j) r,^ 

Mr: » = » = 1) = ^(-1)* E (-D'^^i • («) 

Any three-dimensional overlap is reduced to the one-dimensional form of Eq. ([II]) in a direct 

way 

((/3;R)|(a;r)>= J[ {{(5 x -R x )\{a x -r x )) . (45) 

x=l,2,3 

Similarly we introduce 

((/3,r)|e- lk - x |(a,r')) = r / 3 / 2 x 
/ —k 2 (r — r') 2 v 2 w 2 —ik(rv 2 + r'w 2 ) 

6XP l^^ + w 2 ) ~ 2(v 2 + w 2 ) {v 2 + w 2 ) 

where 

^a(( r '- r )' k ;^ w ) = II ^a.^-r^^^w) . (47) 

x=l,2,3 



^ 5 ((r'-r),k; 



(46) 



B. Calculational details of the A + A — > 2A + tt + reaction. 



Overlaps of many-body nucleon wave functions can be expressed in our approximation 
by a determinant of single-particle overlaps: 



,&;R)|(5 1 ,...;r)(...;r')> 



((A;R)|(ai;r)) 



<(/3 n ;R)|(a i; r)) 



<(/? i; R)|(a n ;r')> 



<(/3 n ;R)|(a n ;r')) 



(48) 



Eq. ( 4q) allows one to take identical exponential multipliers in each row outside the 
determinant as a common factor in all calculations leaving only the matrix of polynomials 
Pz a to be evaluated. A simple example of this is the calculation of the normalization: 

((ax, a 2 . . . a A - v) v ) \(a u a 2 ... a A ; r ')„> = | |P(r -r';v,v)\ \ e -Mr-r')^/^ e < P (r-r') (49) 

In this expression ||P(r — r'; is a determinant of a matrix with the entries Pdidy As 

discussed in Sect ||, this overlap is equal to that of the CM wave functions of two harmonic 
oscillators located at r and r'. For a nucleon system in the lowest state (in terms of harmonic 
oscillator shell excitations), the CM wave function is the harmonic oscillator wave function 
of the ground state |(0,r) f ). We obtain an interesting mathematical fact 



\P{r-r';v,v)\\ = P 00 {r -*';<;,<;) = 1. 



(50) 



Comparison of the exponents in Eq. (|49|) and Eq. ([46]) gives the value of the oscillator 
parameter for the center-of-mass oscillation as q = v\/A. 
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With the same strategy, one can approach the calculation of the reaction A+A — > 2A+tt + 
extracting all exponential factors. Corresponding values of the overlaps F n ,Gi and H n i may 
be rewritten, defining new polynomials /„, gi and h n f. 



Fn = rf(2Z-i)/ 2 x 



exp (-r]vw [Z(R - r) 2 + (Z - 1)(R - r') 2 ] /4) /„ n < Z , 
exp (-r]vw \(Z - 1)(R - r) 2 + Z(R - r') 2 ] /4) /„ n > Z ; 



(51) 



-ffni = ^ 3/2 exp 



= ?f( 27V )/ 2 exp 

- r) 2 v 2 w 2 + k 2 + 2ik ■ (Rv 2 + rw 2 

2(v 2 + w 2 ) 



(R-r) 
\ 



(R-r 



f\2 



/4 91 



(if n > Z, r <^ r') 



It is useful to notice here that all the polynomials are functions of distances between 
the nuclei (r — R) and (r' — R) that we will denote as x and y respectively. Considering 
integration in Eq. (|TBJ) over variables x, y and R we observe from Eq. ( HTJ ) that the 
oscillating phase has the form 



exp 



— ik ■ (Rv 2 + rw 2 ) 



+ zp(r — r') — ipj ■ R 



(v 2 + w 2 ) 

( ik ■ x.w 2 . \ 

= exp hp- x-y — — - - i(k + pf) ■ R , 

y v z + w z J 

and integration over R gives a momentum preserving 5-function that requires k = — P/. For 
convenience we split the sum in Eq. ([16]) over n < Z and Z < n < 2Z and substitute F, G 
and H from Eq. flBTD 



(F\H\I) 



1 



V 



AfiAff J \J2ojV 2m 



^3A e -Avwr,(x 2 +y 2 )/A e ~k 2 /2(v 2 +w 2 ) 



-iV.--x.rjw/2v 



J2 JiQiK + e 



-ik-y-qw /2v 



E (-l) i+j fi9Ai e-^-^dVy . (52) 



i<Z,j 



i>Z,j 



The terms fiQjhij are again some polynomials of x and y proportional to \k\ and containing 
parameters v and w. The final integration can be performed with the help of Eq. 
corresponding parameters a and b being 



a = Arjvw/4 , b = ±p — ki]w/2v . 
As a result, we arrive at the formula fllTD with polynomials 



(53) 



P{k,p) = 
Q(k, P ) 



1 

Ikl 



1 

Ikl 



i<Z,j 

E (-ly^figjhj 

i>Z,j 



— I 



p + kr/w /2v ip \ 
y/Arjvw ' y/Ar/vw I 



( ~ip ■ P — ki]w/2v s 
\ V 'Arjvw ' y/Ar]vw 



(54) 



where the first argument is the transformation of elements of vector x and the second that 
of vector y. From here it is also seen that if before transformation there existed a symmetry 
between x and y, i.e. the nuclei were in an identical state, then P(k, — p) = Q(k, p). 
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C. Toward a complete analytical answer, reaction A + A — > 2 A + 7r°. 



As it was pointed out in the main text, the amplitude of the pionic process is approxi- 
mately proportional to the amplitude of the fusion reaction. One can study the properties of 
the determinants arising in a fusion reaction in a quite general way, separately considering 
the four types of particles distinguished by spin and isospin in the reaction of fusion of the 
type A + A — > 2A. This leads to the following form of a single-particle overlap matrix 



f /<(&;R)|(a i; r)> 



2A I 



<(/3 i; R)|(dU;r')>\ 



\((A4;R)|(« i; r)) (0 A ;R)\(a A ;r')) 



(55) 



first A nucleons second A nucleons 

Without loss of generality, R can be set to zero. A second important feature is that in 
nuclei under consideration all inner shells are filled. Therefore, the resulting determinant 
is a function of a nucleon number A and extra parameters arising from different ways to 
distribute the particles in the outer shells. 

It is interesting to present the exact result for the one-dimensional case where the problem 
is uniquely defined. We consider two oscillators with single-particle states from till A — I 
overlapping with one larger oscillator with occupied states from up to 2 A — 1, 



((0;0) v \(0;x) w ) 

((2A-l;0) v \(0;x ) w ) ■ 
{-l) A {w{x-y)) A2 ^2Ay. 

2^^(11^) 



((0;0) v \(A-l;y) w ) 



{(2A-l;0) v \(A-l;y) t 



7] 



(2A-1)A _ \, 

7 



W(X 



y)) q v~ 



(56) 



The result is just a single term which depends only on the distance between the two initial 
oscillator locations raised to the power equal to the difference in total number of quanta 
between initial and final systems, q = A 2 . The term 77 = 2vw/(v 2 + w 2 ) comes in the power 
of total number of quanta in the final nucleus, Q = (2 A — I) A. This remains true only 
for Fermi systems in the ground state, i.e. if there are no gaps in the harmonic oscillator 
single-particle level occupation. The situation for a three-dimensional oscillator is similar. 
The required polynomial is still given by one term that has a form of the product 

^(x x - y x ) qx (x y - y y ) qy (x z - y z y'w q * +q » +q ' v Qf , (58) 

where integers q x , q y and q z are differences of the number of quanta between the final and 
initial systems in x, y and z directions, respectively. A specific three-dimensional complica- 
tion arises from the following aspect. The lowest energy state is, in general, degenerate as 
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for non-magic nuclei one has the freedom of placing several particles into (n + l)(n + 2)/2 
degenerate levels of the n-th shell. The numerical parameter 7 depends in this case also on 
the way the particles are placed in the outer shell of each nucleus. The harmonic oscillator 
symmetries in the problem often prohibit the transition. 

The polynomials in Eq. (|5"2"D acquire a form of a product of four components, each of the 
form of Eq. (|55D for each type of nucleons, times the sum of terms (a ■ k) acting on every 
pair of interacting nucleon species. Using the integrals from Eq. Q55|) and writing the action 
of (ex • k) between initial and final spin parts of the wave function as a matrix element M 
we arrive at the expression for the polynomial in Eq. (^2]) 



^-l)!!( ft ,-l)!!2*T„(i P ^^j)M . (59) 
In the above expression we have redefined 7 as a product of 7's for all four types of nucleons. 
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